On the Integrability 
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In this letter we present an analytic evidence of the non-integrability of the discrete nonlinear 
Schrodinger equation, a well-known discrete evolution equation which has been obtained in various 
contexts of physics and biology. We use a reductive perturbation technique to show an obstruction 
to its integrability. 
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INTRODUCTION 

The nonlinear Schrodinger (NLS) equation 
i9J + a,,/ = a|/|V, f = f{x,t), a = ±1,(1) 

is a universal nonlinear integrable partial differential 
equation (PDE) for models with weak nonlinear effects 
[l| . It has been central for almost fourty years in a large 
variety of areas in sciences and it appears in many phys- 
ical contexts, see for instance 0, S B S @|- From the 
integrability of this PDE it follows the existence of in- 
finitely many symmetries and conservation laws, and the 
possibility of solving its Cauchy problem, once the initial 
data are prescribed. In correspondence with its symme- 
tries one finds an infinite number of exact solutions, the 
solitons, which, up to a phase, emerge unperturbed from 
the interaction among themselves. 

Many physical and biological applications involve lat- 
tice systems. In literature, one may find a few discrete 
forms of the NLS equation. The most relevant one is the 
discrete NLS (DNLS) equation 



idtfn 



2/n + /n-l) =^|/„|'/„ 



(2) 



where n e Z is a discrete index, t is a, real variable, /„ 
is a complex function and ft, is a real parameter related 
to the space-discretization. Its continuous limit {h —^ 0, 
ri — > oo, x = nh finite) goes into the integrable NLS equa- 
tion ID). Eq. ([2]) is one of the most studied lattice mod- 
els (see for instance 0, 0, 11, 12| and references 
therein). Among the many recent applications of Eq. 

let us just mention the theory of Bose-Einstein con- 
densates in optical lattices [l2| and semiconductors 11 1 . 
The DNLS equation possesses exact discrete breathers 
solutions, where the bumps are spatially localized and 
periodic in time l3, l3|- However, just a few number of 



schemes have been used to exhibit its chaotic behavior 
[l5| . As far as we know no proof of its non-integrability 
is known and a few articles can be found on this subject 

An integrable discretization of the NLS equation has 
been found by Ablowitz and Ladik [l3|. It reads 



ic't/n + - 2./ri + fn-l) 



\fnfifn+l + fn-l) ■ 



(3) 



conserved quantities is known and thus the DNLS equa- 
tion ^ is supposed to be non-integrable [l4| . Numerical 



Eq. ([3]) has an infinite number of explicit standing as 
well as travelling soliton solutions [1]. In the limit as 
h —> 0, also Eq. ^ goes into the NLS equation ([1]). Eqs. 

and ^ exhibit very different responses to the same 
initial data. In the case of an infinite lattice with rapidly 
decaying boundary conditions. Eg. ([3]) has soliton solu- 
tions while Eq. (HD does not [ll, ^ . 

Multiscale perturbation techniques [ll] have proved to 
be important tools for finding approximate solutions to 
many physical problems by reducing a given PDE to 
a simpler equation, which can be integrable [l[. These 
multiscale expansions are structurally strong and can be 
applied to both integrable and non-integrable systems. 
Zakharov and Kuznetsov [20| have shown that, starting 
from an integrable PDE and performing a proper multi- 
scale expansion, one may obtain other integrable systems. 
In particular, they showed that the slow-varying ampli- 
tudes of a dispersive wave solution of Eq. ([T]) satisfies 
the Korteweg-de Vries (KdV) equation ([T]) and viceversa. 
Calogero et al. [2lj have used the multiscale perturba- 
tion technique as a tool to give necessary conditions for 
the integrability of large classes of PDE's both in 1+1 
and 2-1-1 dimensions. In particular the non-integrability 
of the resulting multiscale reduction is a consequence of 
the non-integrability of the ancestor system. Multiscale 
techniques have been used also in the papers 22, 3, 24 1 
to prove integrability of several PDE's. 
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Some attempts to extend this approach to discrete 
equations have been proposed 

32 , 3^ . In S^, 31 , 3^ one can find a multiscale expansion 
technique on the lattice which, starting from dispersive 
integrable Z^-lattice equations, provides other Z^-lattice 
equations. To do so one had to introduce a slow-varying 
condition on the amplitudes by requiring 



(A„)^+Vn = 



(4) 



p being a positive integer and A„/„ = /„+i — /„. As a 
consequence, the resulting reduced equation turned out 
to be non-integrable even if the ancestor equation was in- 
tegrable. However, as shown in [27|, Hp = oo the reduced 
equations become formally continuous and their integra- 
bility may be preserved by the discrete reductive pertur- 
bation procedure. In this way the multiscale expansions 
easily fit with both difference-difference and differential- 
difference equations. The illustrative example considered 
in [i^l has been the lattice potential KdV equation, a dis- 
persive nonlinear Z^-lattice equation obtained from the 
superposition formula for the soliton solutions of the KdV 
equation. There one performed the multiscale expansion 
of the weak plane wave solutions of the discrete dispersive 
linear system. A proper representation of the discrete 
shift operators in terms of differential ones, provides the 
integrable NLS equation ([T]) as the lowest order sccu- 
larity condition from the multiscale expansion. Further 
examples have been considered in (29| . They confirmed a 
discrete analog of the Zakharov-Kuznetsov's claim [iot : 
"if a nonlinear dispersive discrete equation is integrable 
then its lowest order multiscale reduction is an integrable 
NLS equation". 

In the present letter we consider the multiscale per- 
turbation analysis of Eq. ([2]). Multiscale analysis of the 
DNLS equation ([2]) has been already considered in [i^t . 
giving a differential-difference system which does not ful- 
fill any integrability criterium. However this result did 
not prove the non-integrability of the DNLS equation, 
as similar results have been obtained in the case of the 
integrable equation ([3]) [l^, [l^l . Here, extending to lat- 
tice equations the approach used in [i^, |22„ 22e i23| and 
computing the higher order terms in the reductive per- 
turbation expansion we are able to provide an analytical 
evidence of the non-integrability of Eq. ([2]). In fact, even 
if its lowest order multiscale reduction is an integrable 
KdV-typc equation, the higher orders reductions exhibit 
non-integrable behaviors. 



DISCRETE MULTISCALE ANALYSIS OF THE 
DNLS EQUATION 

By the position = [z/„(t)]i/2exp[i0„(t)] the DNLS 
equation ^ may be written as the following system of 



real differential-difference equations: 
c^tt'n + ^ (^VaI"sin/3+ -f Vojf sin/?;^^ = , (5) 

^ , /tt^cos/?,;^ + (6) 

0, 



where a^{t) ^ i^n^-nii, P^(t) = (f>n±i - 0„, and 7^(0 = 
i^n^^n±i- By analogy with the continuous case, see [20|, 
we expand the real fields Vn{t) and (jjnit) around the con- 
stant solution fn{t) = exp {—lat) of the DNLS: 

oo 

i^nit) = 1 + ^ e^^ {t„,}„,>i) , (7) 

oo 

0„(O = -at + Y,^^'''<l>^'H^At,n},n>l), (8) 

2=1 

where e, with < |e| <C 1, is the perturbation parameter. 
The fields iy^^'> and in Eqs. dUE]) depend on the 
slow-space variable k — eQn, ^ S M, and the slow-time 
variables tm = e^™~^i, m > 1. The free parameter C, can 
be fixed so as to obtain the proper continuous limit. 

Given a function = v{k, {tm\m>i) we expand 

Un±i{t) and dtUn{t) in terms of the slow variables k and 
>i. Introducing the shift operator r„ such that 
T^Un = Un±i we have 



Un±l{t) = (7^) v{K,{tm}m>l) = 

OO 

= ^ eMf {<m}m>l) : 



i=0 



dtUn{t) = ^e^' ^dt,v{K,{t„i},n>l) : 



(9) 

(10) 



where, defining = (T^ — 1)* to be the i-th order dif- 
ference operator, one get 



^ 7 



(11) 



i=i 



If u„ is a slow- varying function of order p, see Eq. Q , 
we can truncate the infinite series in Eq. (jlip . In such 
a case the ^K-operators reduce to polynomials in the A^- 
operators of order at most p. Hereafter we shall assume 
p = oo and the ^^-operators are formal differential oper- 
ators. 

Taking into account the expansions ([7][S|) and Eq. © 
we have the following formulas for the shifts of the func- 
tions Vnit) and 0„(i): 

+ 00 b72] 

'^n±l{t) = 1 + E E ^'Af_,,iy^'HK, {tr,U>l) , (12) 
j=2 i=l 

4'n±i{t) = -a-t + 

+00 [b + l)/2] 

+ E E ^'Af_,,^,cj^^'H^dtrr.}rn>l), (13) 
j=l i=l 
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where [x] denotes the integer part of x. 

Let us introduce the multiscale expansions (|12m3p . to- 
gether with Eq. pU]) . into Eqs. ([5][n]) and require that 
these equations be satified at all orders in e. 

At the lowest non-trivial order one finds i^*^^-' = 
—(jdt-i^(j)^^\ From now on all results will be presented 
just for the functions 4>^^\ 

At order we get (d^.^ — c^5^ 4>''^^ = 0, where c = 
± (Ccr^/^) /h. As c has to be real then (7 = 1; moreover 
we choose C, = h so that c remains finite as ft, — > 0. 
Therefore the resulting equation at this order is satisfied 

by (/.(I) = {tm}m>2) with e = K - C<1. 

At order e^, the no-secular term condition implies 
- c^J^) 0(2) ^ 0, i.e. ^ 0(2)(^,{t„,}„>2). The 
evolution equation for 0^^^ w.r.t. the slow-time t2 reads 



K2 



(14) 



a|0(i)-l(.,0(i))^ 



with a = c(3 - /i2)/24. Eq. ^ is a potential KdV 
equation. Therefore, if Eq. ^ has to be integrable, then 
its multiscale reductions should provide the integrable 
evolution equations (j > 3): 



C[fm=dif{S,) 



920(1) 



(15) 



2a 



dyf{y) 



where C is the recursive operator associated with the 
KdV hierarchy and the 6j's are free coefficients to be 
fixed later. 

We assign a formal degree to the ^-derivatives of the 
fimctions 0(^), deg(5|0(j)) = I + 2j - I, t > 0, and 
define Vn as the vector space spanned by the products 
with total degree n of all derivatives 9|0(''). We denote 

hyVn C Vn the subspace spanned by those products of 
derivatives 9|0(^) with j < r. Similar vector spaces have 
been introduced in 2^, 2^ in the case of the multiscale 
analysis of PDE's. 

After caring for secularities, the order yields 0(^) = 
"Z*'^) (Ci {^m}rn>2) and thc following non-homogeneous 
evolution equation for thc field 0(2) w.r.t. the slow-time 
^2, depending on 0(i) and its derivatives: 

dt.^^'^ - ^i^^9|0(2) + ^a,0(i)a,0(2) = (16) 



5/i2-7 



(a|0(i))V^(a,0(i))^ 



64 V"^^ 



16 ^ « 1920 



Substituting Eq. into Eq. HH) with j = 3 and 

fixing 63 = —c{h^ — 30ft2 — 15)/1920 in order to remove 



residual secularities, Eq. (flBl) reduces to the following 
evolution equation for thc field 0(2) w.r.t. the slow-time 
t2: 



dt. 



K'2 



,(2) ^ J{t2) 



(17) 



where K'^ (0(i)) 0(2) = (d/d6')X2(0(i) + 6(j}(^^e=o is the 
Frechet derivative of K2 (0(i)) along the direction of 0(2) . 
In Eq. (|17p the forcing term /(*2) is a well-defined ele- 
ment of 'P'^\ namely a linear combination of three in- 
dependent differential monomials, with coefficients that 
are rational functions of h. At this same order, the inte- 
grability of Eq. ^ implies the existence of the following 
evolution equation for thc field 0(2) w.r.t. the slow-time 



dt 



,0(2)-/^^(0(l))0(2)=/(*3) 



(18) 



where /(*^) is an element of the space v'^\ which 
if Eq. ([2]) has to be integrable, must satisfy 
the compatibility condition [St., — K'.^ {4>^^^y\ /'•*^) = 
[dt^ - K'^ (0(1))] /(*3). Such a condition allows one to 
express the coefficients of the polynomial /(*^) in terms 
of those of /(*^), and does not impose any further con- 
straint on the coefficients of /(*^). As this condition is 
satisfied, eventual obstructions to the integrability of Eq. 
([2]) will appear at higher perturbative orders. 

Let us now consider the order e^. The resulting equa- 
tions provide the evolution of the field 0('^) w.r.t. the 
slow-time t2. This is given by an integro-differential 
equation. Introducing the fields (/3(^) = 9^0(-'), taking 
care of secularities and taking into account that 0(i) 
evolves w.r.t. the slow-time t4 according to Eq. (jlSp . 
we get 0^^) = 0('i)(C, {tm}m>2) and 



2 [f 



0(1; 



(3) 



,(t2) 



(19) 



where H2 (v''^'') v'^' is the Frechet derivative along ip'-^'^ 

of thc KdV flow H2 (<^(i)) = d^K2 (0(i)). Here g(*^) e 
(2) 

T'g is a linear combination of fourteen differential mono- 
mials, whose coefficients are well-defined rational func- 
tions of h. The evolution equation of (^(^) w.r.t. the 
slow-time ts takes the form 



(3) 



<^(i' I p 



.(3) 



9 



(t3) 



(20) 



where i^^^^) f^^^ is the Frechet derivative along tp^^'> 
of the higher order KdV flow H3 ((^(1)) = 8^X3 (0(i)). 

(3) 

Here gl^ is an element of the 31-dimensional vector space 

(2) 

Vli whose coefficients are determined by requiring the 
compatibility condition 



o(l) 



i?^ ^(1) 



Ms) 



(21) 
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Eq. (|2T|) is a necessary condition for the integrability of 
Eq. ([2]) . In this case only nine out of the fourteen coef- 
ficients of 9'*^' are independent. Thus we have five in- 
tegrabihty conditions, whose explicit expressions will be 
published elsewhere in a more detailed paper [3J|. The 
obtained constraints on the polynomial g^*^^ are not sat- 
isfied by the coefficients computed in Eq. ([19]). Then, 
this incompatibility implies that the DNLS equation ^ 
cannot be integrable. 



CONCLUDING REMARKS 

By performing a discrete multiscalc analysis, wc have 
proven that the DNLS equation ([2]) is not integrable. Al- 
though its lowest order reduction gives a KdV equation, 
the higher order terms do not satisfy the required integra- 
bility conditions. It is remarkable to note that a similar 
analysis performed on the integrable discrete NLS equa- 
tion ([3|) provides integrable reductions [sl ]. 

Moreover, by using the expansions (I7][5]), our pertur- 
bation technique enables us to construct approximate so- 
lutions of Eq. ([2|) starting from the exact solutions of the 
integrable Eq. (11411 . However these solutions will break 
down in the far-field region, due to non-integrability of 
Eq. (HD. 
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